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ABSTRACT . , . . 4. ^^n 

The VARAN (variance Analysis) program is an addition 
to a series of computer programs for, multivariate analysis of 
variance. The development of VARAN exploits the full linear model. 
Analysis of variance, univariate and multivariate, is the program s 
main target. Correlation analysis of all types is available with 
printout in the t-ernacular of .correlation. The hybrid of these, 
homogeneity of regression, has been added with as much flexibility as 
can be currently mustered. In addition to these, VARAN includes 
several styles of factor and component analysis complete with tests 
of factorization and rotation techniques. This research memoranduB is 
the manual for the second' editio^ of VARAH, an enlargement of the 
first edition. The mainstream of the program is essentially unchanged 
but several additions have be^n made and three small programming 
errors have been corrected. The most extensive addition has been in 
serial correlation analysis. (Author/RC) 
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Introduction 

< 

'mis memorandum is the manual for the Vmm program. (The acronym 
comes from the words VAHiance ANalysis) 

The VARAN program is the latest" addition to a series of computer 
programs lor multivariate analysis of varia^ .e which originated a'bout 1957- 
About that time R. Darrell bock had a MANOVA type program running on the 
UNIVAC at the University of North Carolina. Bock's original program 
handled only a few variables and was otherwise quite restricted and was 
not in much general use. Later, in I962, C. E. Hall and Elliot Cramer 
with Bock's assistance published a program called MANOVA in FORTRAN 2 which 
was internationally circulated. The program handled 25 variables and 100 
degrees of freedom for hypothesis. 

The wide use of this program prompted further development by Cramer 
and by Bock and Finn. In I966 Cramer published a program called MANOVA 
in FORTRAN if followed by Finn and Bock's program MULTIVARIANCE also in 
FOETRAJi 4. The development of these two programs greatly increased the 
scope of multivariate analyses which could be performed on computers. 

NRJLTIVARIANCE was the first of the programs in this series that 
utilized the full linear model. The earlier programs had been restricted 
to models of the cell means in their main streams of calculation. Other 
variations on the linear model, like canonical correlation, were of an 
"accidental" nature. With the development, of MULTIVARIANCE, the series 
turned to the flexibility of the full linear model. 

The development of VARAl^' continues this series in the exploitation 
01 the full linear model. Analysis of variance, univariate and multivariate, 
is the main target of the program as with the earlier programs. Corre- 
lation analysis of all types is av&i].able with printout in the vernacular 
of correlation. The hybrid of ohese, homogeneity of regression, has been 
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added with as much flexibility as can be currently mustered. In addition 
to these, VARAN includes several styles of factor and component analysis 
complete with tests of factorization and rotation techniques. 

The addition of factor analysis a.nd correlation techniques to a 
MANWA program brings several new capacities to multivariate analysis in 

f 

addition to the customary univariat . a^d multivariate ANOVA, correlation 
and homogeneity of regression problems. 3ome of the new capacities for 
analyses are listed below. 

1. Three kinds of batt.ery red^ ?tion procedures from multiple 
correlation applied to discriminant analysis, multivariate analysis C'f 
variance and canonical correlation. 

2. Seven kinds of r9t?.tion schcfn^s from factor analysis applied to- 
discriminant analysis, multivariate analysis of varianc- and canonical 
correlation . 

3. Factor extension from factor analysis applied to discriminant 
analysis, multivariate analysis of variance and canonical correlation. 

4. Estimation of population varia-ce and covariance parameters 
from several samples as I'ou^d ]n analy^'. of variance applied to corre- 
lation and factor analysis models. 

5 Use of "dummy parameters^' from analysis of variance applied to 
factor analysis, analysis of /arianjc (obtaihinc intraclass correlation 
coefficients) and correlacion analysis (complex biserial correlation). 

6. Homogeneity of regression te^-hniiaes from analysis of covariance 
applied to correlation analysis from multiple samples and factor analysis 
from multiple examples. 

Y. Iriteraction tables for ANOVA models. 

8. Variance reduction analysis for determinin/:^ the effects of non- 
orthogonality in ANOVA analysis. 
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9. Dimension reduction analysis for determining overlap of effects 
in MANOVA and complex canonical correlation analyses. 

The VARAN program was written more in the style of a mathem'^tical 
exercise in linear calculations than as a solution for specific statistical 
models. Therefore, the user can expect to find applications in analysis 
which were not specifically contemplated by the authors. The VARAN j)rpgram 
was also constructed to be easily expanded to include new linear model 
techniques as they are generated. Updates will be forthcoming periodically 
and for this reason it is suggested that all copies be obtained either 
from Educational Testing Service or from the authors. 

As is customary with programs of this size, the. author makes no 
ironclad claims of arithmetic accuracy. Any errors discovered by the 
users will be quickly corrected and distributed to other users. It should 
be noted that there were 85 or more problems run to check the accuracy 
of the main streams of calculation. It may also be noticed that the 
ability to make linear transformations of the data provides a wide variety 
of internal checks on calculation accuracy. 

All arithmetic is single precision except for orthogonal polynomial 
construction. The user is warned thp.t lengthy manipulations of highly 
correlated variables are not advised. (This is not much of a ieficiency 
in handling ANOVA designs since data are not generally useful when highly 
correlated and cell counts are generally close to orthogonality.) 

The authors wish to express their appreciation to the programming 
staff of Educational Testing Service for three years of assistance. 
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int reduction to the Second Edition 

The second edition of VAEAN is an enlargement of the first edition. 
The main stream of the program is essentially unchanged but several adr'itions 
have been made and three small programming errors have been corrected. (No 
errors have been found in the arithmetic.) 

The most extensive addition has been in serial correlation analysis. 
I'his addition should prove useful in sociological surveys and in complex 
sampling designs. 

A feature has been added which makes Potthoff and Roy models easier 
to use. 

An orthogonal Procrustes rotation scheme has been added to aid specific 
hypothesis testing in canonical analysis. 

Descriptions of three machine-dependent subroutines have been added to 
the manual to ease conversion to other machines. 

The test problems have been reworked to provide tutelage as well as 
checks at new installations. 

The manual has been added to and revised only slightly for VARAN2. 
None of the program set up has been deleted although additions have been 
made. The old manual still works for the new program. 

It is hoped by the authors tnat these additions will make VARAN2 even 

more useful. 



7 

ERIC 



ERIC 



Control Cards in Brief 

1, TITLL Cards : (optional) 

Cols 

1-4 TITL 

>-80 Any alphameric description 

As many cards may be used as desired 

f 

2. PROBLEM Card : (required) 

• Cols 

1-U PROB 

5-6 Number of variable format cards for data (lO or fever) 

8 Number of contrast card sets (at least 2, at most 7) 

10 Number of individual '^significance test control cards 

(at most 6) _ / 

12 ' Print cell means and variances (1 means yes)* 

Ik Print reduced model matrix (l means yes) 

/ 

17-18 Data file number if not cards 

20 ' Print only estimates or rav regression coefficients 

(1 means yes) 

22 Controls printed output (0 gives minimal output). 

See Table I 

2k Controls printed output (0 adds nothing to above), 

^ee Table I 

26 Type of covariance adjustment (0 or 1 does the 

classical adjustment: 2 -ises "error" regression 
weights for reduction) 

28 ' Do not copy data input file for reanalysis 

(l means don't) 

3* Contrast Card Set : (At least tv;o required) 
Cols 

1 Letter identification of the factor (W is not allowed); 

V is allowed only for continuous variables 

Q blank always means "no'' 

B 
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Cola 



^-3 



k-6 

8 

10 



2 col 
fields 

Next 2 col 
field 

Next 2 col 
field 



Next 2 col 
field 

20 col 

field 



Contrast code for design variabler^ 

= 0 Regular contrasts (also called nominal or 
deviation) 

- 1 Special contrasts 

= 2 Orthogonal polynomials (integral values) 
(19 levels or less) 

= 3 Reverse Helmert contrasts 

= k Special design parameters 

= 5 Orthogonal polynomials (normalized values) 
(39 levels or less) 

Function code for continuous variables (after V) 

- 0 do nothing 

= N obtain Nth power of all variables (N > l) 

Number of levels of a factor (at most kO) or, after V, 
the number of continuous variables (at most 39) 

Read in variable names (l means yes) 

Controls linear transformation of the variables 

' = 1 Linear transformations of the variables 
to be read in 

= 2 Weighting matrix for transformations 

followed by linear transformations (Potthoff 
and Roy models) 

= 3 Inverse of weighting matrix for transformations 
(Potthoff and Roy models) 

Number of variables in each partition (if no partitioning, 
do notning) (at most NORD partitions) 

(two commas required) 



Recoding of ail tht. factor identification or, after V, 
reordering of all the vari£.jles (if no recodirig, or 
reordering, do nothing) 

(two commas required) 



Any description of the contrast or variables 



1 col 



. (a period required) 
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k. Contrast Card Set Options: (optional) 



a. Variable names cards. FOW'IAT (10A8) 

Any names in 8 column fields, 10 names per card. The name cards 
must precede other options . 

b. Special contrasts matrices or special design parameters. FORMAT 

(Bfio.o) 

Contrast or design matrices are entered a row at a time, the row 
for the grand mean being first. Contrast and design matrices are 
square and of dimension equal to the number of levels of the lactor 

c. Orthogonal polynomial metric. FORI-iAT (8FIO.O) 

Ihe coefficients of the linear polynomial fit only. 

d. Linear transformations and weighting matrices. FORMAT 
(8F10.0) 

Weighting matriced for transformations (col 10 = 2 or 3) precede 
the trans formatioi^s. Each row of a weighting matrix is pionched 
in 10 column fielfs, 8 entries per card on as many cards as 
necessary. Each -row must begin on a new card. The weighting 
- matrix must be nonsingular, square and of the order of the number 
of variables . 

For linear transformations (col 10 = 1, 2 or 3), each regression 
coefficient is punched in a 10-column field, 8 coefficients per 
card, as many cards as necessaxy for each regression equation. ) 
Each equation begins on a new card. The transformation matrix 
must be nonsingular square and of order equal to the total number 
of continuous variables. 

5. Significance Test Card (required) 

This card contains the model statement and indicates the tests io be 
calculated. T^e model statement is limited to at most 9 cards. Some of 



the more common models are these. 



a. Two way factorial analysis of variance 
ANOV A : MS , A=V , B=V , AB=V . 

b. Correlation analysis 
CORREL:\^0,V1=V2. 

c. One way analysis of covariance 
ANOVA:Vf=0,A;=Vl/V2. 



d. Minres factor analysis with varimax rotation 
FACTOR^: VMS, 0=V. 

e. Interaction Tables in a three way factorial 
INTA:BL : WAB=V , WAC=V , WBC=V . 

f . Serial correlations between two categories, B, pooled 
across two samples, A. 

SERCOR : W=0 , A=0 ; BWA=V . 

The following symbols are used: 

ANOVA, CORREL, FACTOR, SERCOR, INTABL and VARED are acronyms 
to begin the card and state the kind of model 

I 

: to separate acronym from model statement 

Letters are used to denote ANOVA-type "main effects" 

Letter folloved by a numeral for partitions (except after W) 

A sequence of letters or letters and numerals without symbols between 
for interaction effects 

-h to indicate pooling of ^effects 

= to separate hypothesis variables from error variables. 

In this (description any variable or variables 
to the left of an equal sign is called a "hypothesis" variable 
any variable or set of variables to the right of the equal 
sign is called an "erroi^' variable except if it is also to the 
right of a slash (/) or an ampersand (&). 

- (minus) for dimension reduction 

^ (comma)"to separate tests' 

; (semicolon) in serial correlation: Any hypothesis variables 
to the left are to be treated as contrasts among sample means, 
any hypothesis variables to the right are to be treated as 
contrasts among categories 

/ to indicate that covariates follow 

& to indicate that extension variables follow 

* in an error term between two sets of variables to 
indicate that the regression of the left set on the 
I right set is to form an error tern, for an analysis 
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0 (a zero) as an error term indicates that the hypothesis 
term is to be included in the model but not tested 

0 (a zero) as a hypothesis terra" IHTilcertes that the 
error term is to be used in a factor analysis 

W as a hypothe.^is te>T. ih'licatt j the cor.Jr'tant term 
or the ^^rand mear \^ 

W ii: a letter-number sequence LridiOaTtes that the effect 
on the left is nested in the effeci: or ^eJls on the 
/right 

' ^ di^it is required arter K'vCTOR to indicate the type 
' of factor analysis, see Table II 

\ after FACTOH and its di^^it indicates that rotation of 
^ Uie factors Is to be done: a digit M^^ust-follovr this to 
indicate type of rotation^ see T-^ble III 

. (a period) must ena the statement 

iC\<JIi'. denotes an analysis of Variance model, multivariate or univariate; 
(JOi<i^j^]L ler.otes a correlation model, product monent, multiple or canonical; 
SERCOH denotes a serial correlation model, biserial or canonical; 
I'ACTO.. de :)tes a far-tor analysis model, including principal conponents; 
l^iTABL denotes Inl-raction tables for analysis of variance; VARED denotes 
a variance reduction study of correlation and analysis of variance problems. 

Variable Format Cards (required) 

Tl-iese cards describe the data records and are the usual variable format 
cards of FORTRAN. Ten cards are allowed. (There is no fixed format data 
option in this program.) It is required that cell identification be 
rea<i ir. oefore the scores on the continuous variables. 
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7, Data (required) 

The data may be on cards or a^ computer kept file (see columns l/'-lB or 
the PROB card). If the data are on cards, put en v s cards after tht 
data to make a blank, but complete, data record. If the data an- on a tape 
or disk file, the end-of-f ile mark will signify the end of the data. 

8. Individml Significance Test Control- Cards (optional) 

These cards control the use of expository calculations on designated 
tests of the model. At most 6 of the significance testk of the model ca- 
be subject to these techniques. 

At 

Cols . Information 



1-2 The number of the significance test to vhich this infor- 

mation applies (obtained by counting equal signs from 
the right [or the period]) 

3-1*. Probability statement "cimes 100 for battery reduction, 

dimension reduction and rotation of canonical variates 
in ANOVA arfd CORREL 
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Cols Information ■ ^ 

6 Type of battery reduction procedure for error variables 

1 for Efroymson's stepwise addition/elimination 
procedure 

2 for Wherry-Doolittle stepwise addition procedure 

3 for Step-up variable elimination procedure ^ 

7_8 Number of factors to be obtained in a factor analysis. 

For components analysis, the number of components is 
always the number of variables and for principal factor 
analysis this is always one less than the number of 
variables; no entries are required for these solutions. 
All other solutions require an entry here. 

IQ Factor' analysis tests 

Blank for none 

1 for Rao's test in canonical factor analysis 

2 for Rippe'-j test 

Type of rotation technique (see Table Hi) 

1 for direct rotation of canonical variates; all factor 
analysis 

2 for indirect rotation of canonical variates 
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2£ 1 for taxonany of variables 

2 for taxonomy of groups 

Tft Power of 0.1 (multiplied by lO) for the convergence 

criterion in uniqueness iterations of canonical factor 
analysis (sec TabK- IV). Entering l6 produces a 
criterion of u.ll.6= .0251189 



1 if squared communaiities are supplied for factor analysis 



V P^,, 1 if wciehts are supplied for rotations of factors or 

canonical variates 

9. Individual r .lf ^nificance Te st Card Options (optional) 

a. Communalities for factor analysis. FORMAT (8F10.0) 

Squared communalities are entered in 10 column fields, 8 per card, 
for as 'many cards as necessary. 

b. Weights for Rotation of factors and canonical variates. FORMAT 

(8F10.0) 
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Weights are entered in 10 column fields. For ei ch canonical variate 
or factor there must be as many weights as there are variables. TVie 
weights for each canonical variate must begin on a new card. Wl*en 
rotating canonical variates in ANOVA or CORREL, there must be as 
many sets of weights as^there are canonical variates; if the rotation 
is a taxonomy of groups, each set of weights will have only as many 
weights as there are groups. For principal components, weights 
must be supplied for all conponents (as many as there are variables). 
For all forms of factor a> xlys is, weights must be supplied for as 
many factors as are indicated in columns 7 and 8 of the individual 
significance test card or, if this is zero, as many as the number 
of variables . 

10/ TITLE Cards (optional) 

TITLE cards may be used here to identify reanalyses. 

11, R eanalysis Card (optional) 
Colg 

1-U ANLY (not ANALY) 

6 1 if a new significance test card is used 

8 number of new contrast card sets 

(must be used for changing the partitioning of the 
continuous variables) 



The following features are the same as on the PROB card and appear in 
the same card columns. 

10 Number of individual significance test cards for this 

analysis 

12 Print cell means and variances (l means yet.) 

lU Print reduced model matrix (l means yes) 

Data file number not necessary 



15. 



Cols 
20 

22 

2k 

26 
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Print only estimates or raw regression coefficients 
(l means yes) 

Controls printed output (0 gives minimal output). 
(See Table I) 

/Controls printed output (0 adds nothing to the above). 
(See Table l) 

Type of covariance adjustment (0 or 1-does the classical 
adjustment, 2 uses "error" weights for reduction) 



/a,l other control cards and their options remain the same as for the 
original analysis. 

12. Several Problems 

Many problems may be submitted with a single run. 

13. FINISH card 
Cols 

1_6 FINISH 
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Table 1 
Printed Output Controls 



Information which is printed when column 22 and column 
2U of the PROB or ANLY card are coded. 



y— 

/ 

Codes 
Col 22 Col 2U 


Cross-index 
for Table 

TV 

Information 


Correlation 


1 ^ >* 


Correlation among the hypothesis variables 


a 




Standard deviation of hypothesis variables 






Cross-correlations between hypothesis and error 
' variables 


b 


2 3 


Raw score weights for regressing hypothesis 
variables onto error variables 


c 


2 3 


"Standard score weights for regressing hypothesis 
variables onto error variables 


d 


1,3 


Standard error about the regression line{s) 




b,l,2 


Correlations among the error variables 


e 




Standa d de^dations of the error variables 




t 


Regression sums of squares when total sums of 
squares of error variables are unities 


f _ 




Statistical Summary 






Dimension Reduction Statistics {Canonical Correlation 
only) 



^-unching a blank '(b), a 1 (one) or a 2 (tvo) forces printing of the information. 
**This designates that col 21^ does not, control printing of this information, 
numbered footnotes can be found a-c the end of the table. 
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Codes 



Col 22 Col 2h 

b,2 
1,2 



1,2 

b,2 
b,2 

b,2 



1,2,3 



Irvformation 

Univariate F statistics for error variables 

Correlations between the variaMes and the 
canonical variates within the error set 

Standard score weights for regressing variables 
onto canonical variates within the error set 

Raw score weights for regressing variables onto 
canonical variates within the error set 

/ 

Weights for regressing raw score variables onto 
error canonical variates with unit sums of 
squares 

/ 

Regression sum^of squares when total sums of ^ 
squares of hypothesis variables are unities 

*- 

Univariate F statistics for hypothesis variables'" 

Correlations between the variables and the 

canonical variates within the nyDOthesis set 

Standard score weights for regressing variables 
onto canonical variates within the hypothesis 
set 

Raw score weights for regressing variables onto 
canonical variates within the hypothesis set 

Weights for regressing raw score hypothesis 

variables onto hypothesis canonical variates 
with unit sums of ijquares 

Print only the infovnation indicated in column 
2h 



Analysis of Variance 



b,2 

2 

\ b,2 



Raw estimates of effects 
Standardized estimat^es of effects' 
Orthogonalized estimates of effects 
Univariate standard errors 



o 
n 
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Information 

Error dispersions reduced to correlations^ 

Hypothesis sums of squares when error sums 
of squares are unities 

Hypothesis s-ums of squares and error sums 
of squares (univariate) 

r>tatisti(;al summary 

Dimension reduction statistic^ 

(discriminant analysis and MAHOVA only) 



Univariate F statistics 



\ 



Correlations between the variables and the 
iianonical variates 

t 

Correlations between the variables and the 
canonicaP^^iates weighted by the square ^ 
roots off the associated canonical variance 

Discriminant function coefficients for 
standard scores 

Discriminant function coefficients for raw 
scores 

Estimates of effects for the canonical variates 

Transformation matrix for obtaining canonical 
contrasts (unre'duced)^ * 



Factor Analysis 



Correlation matrix or covariance matrix (as 
appropriate ) 

10 

Estimates of standard deviations 

Comraunality estimates, when applicable 

Final communal itv estimates (canonical factor 
analysis only) 

Characteristic roots and vectors 



IS 
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Col 22 Col 2k Information 

Factor coefficients (unrotated) 

Images and anti images (image analysis only) 

Canonical correlations (Canonical factor analysis 

only) 

I Regression weights for factor scores (total 

. variance procedures only) 

Statistical tests (if usod) 

I Residuals from the correlation matrix 

1 Mean and standard deviation of residuals (below 

the diagonal) 

Frequency /distribution of residuals (if program 

is set, up to handle 20 or more continuous 
variables ) 

Variance Reduction 

^ Sum of squares among design' parameters and 

percentage loss 

* Correlation among design parameters before 

reduction 

Coj^j^elation among design parameters after redaction 

Sujn of squares among contrasts and percentage 
loss 

Correlations aniong contrasts before reduction 

' Correlations among contrasts after reduction 

Correlations among hypothesis sums o^f squares 
before reduction 

Correlations among hypothesis sums of squares 
after reduction 

Data on trace of hypothesis sum-of- squares matrix 

Hypothesis variance data for individual variables 
and percentage loss 



1 
1 

1 
1 

1 
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Codes 

Col 22 Col 2U 



Cross index 
for Table 

Information ' 



Serial Correlation 



2 1 Category variables cross products matrix 

Category frequencies for nominal contrasts ^/^ - 

1^2 Coi^relations among category variables a 

1^2 1 Standard deviation of category variables 

1^2 Cross correlation between category and 

c error variables 
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0,1,2 Estimates of contrasts ajnong category means 

(pooled across samples'*) 

2^ 2 Standardized estimates of contrasts among 

category means (pooled across samples) n 



0,2 Correlati( ns among the error variables ^ 

0,2 Standard deviations of the error variables 

1^2 Regression sms-of-^ueres when total s'oius- 

^ of-squares of error variables are unities a 

Statistical summary 

Dimension reduction statistics 

0^2 Univariate F statistics for error variables g 

Q 2 Correlation between the variates and the 

' canonical variates within the error set h 

2 3 Standard score weights for regressing variables 

onto canonical variates within the error set i 

2 3 Raw score weights for regressing variables 

onto canonical variates within the error 
set , \ ^ 



/ 



Codes 
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Col 22 Col 2k 



0,1,2 
1,2 

0,2 
0,2 



3 
3 



Information 

Weights for regressing raw score error 
variables onto error canonical variates 
with unit svirns of squares 

Estimates of effects for canonical 
contrasts among categories 

Regression sums-of-squares when total 

sums-of-squares of category variables 
' are unities 

Univariate F statistics for category 
variables 

Correlations between the variates and the 
canonical variates within the category 
set 

Standard score weights for regressing 
within the category set 

Raw score weights for regressing variates 
onto canonical variables within the 
category set (transformation to 
canonical contrasts) 

Weights for regressing category variables 
onto canonical variates with unit sums- 
of-squares 
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Footnotes for Table I 



Regression sums of squares when the total sums of squares of the error 
(hypothesis) variables are unities. This matrix has on its diagonal 
the squared multiple correlations between the error (hypothesis; 
variables and all hypothesis (error) variables. The orr-d^agojisl 
elements become the correlations among the regressed variables when 
they are divided by the square roots of the. diagonal elements. 



Univariate F statistics for error (hypothesis) variables-. These are 
the F tests of the multiple correlations for each error (hypothesis j 
variable regressed on all the hypothesis (error) variabDes. 

Standardized estimates are the raw estimates divided by their standard 
errors. A standardized estimate of +1.0 is one standard error above 
the grand mean of the data. These figures are easy to relate to 
confidence intervals abo- the grand mean. 

Orthogonalized estimates are appropriate in nonorthogonal designs. 
They are what is left of the raw estimates after the nonorthogonality 
ff ^he design has been accounted for. The analysis which - P^o--^ 
is an analysis of these estimates. Comparison of the raw estimates , 
with the orthogonalized estimates .s sometimes ^-^^^K'^^lll^Y'''^ \ 
the effects of nonorthogonality on the analysis. In E. Oramer 
MANOVA these are the "Estimates." 

Error dispersions reduced to correlations. In a multivariate analysis 
of variance the error term is a variance-covariance matrix. This is 
that error term reduced to a correlation matrix. 

Hypothesis sums of squares when err6r sums of ^^^3 
tS diagonals of this matrix when multiplied and divided the degrees 
of freedom give the univariate F-ratios. The off-diagonal entries are 
° type of "covariance F" and reflect the relationships among treatment 
effects on the variables, 

correlations between the variables and the canonical ^^^^^^J^^^f ^^'^^ 
by the square roots of the associated canonical variance. These values 
are related to "Student's" t: when the analysis is of two samples and 
a single variable, it is "Student's" t. In a MANOVA or discriminant 
an^Jfis, Bumming'the squares of these values for one variable across 
?Se canonical variates will produce the univariate F ratio for that 
variable . 

F-timate of effects for the canonical variates. These are the mean 
■ discriminant ^c^rcs when the grand mean is zero. These estimates always 
add to z^ro for each canonical variate (or discriminant variable). 
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Transformation matrix for obtaining canonical contrasts. This matrix 
when used to multiply the contrasts, produces the canonical contrasts 
In orthogonal designs it produces exactly the canonical contrasts; m 
nonorthogonal designs it produces the canonical contrasts ignoring 
adjustments for the lack of orthogonality. 

Estimates of standard deviations. This program always assumes that 
the data area sample and not a population. These are not standard 
deviations but sample estimates. The estimates are residual to any 
covarlates or sampling design. 
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Table II 
Factor Analysis Codes 



Codes 


Type of Factor Analysis 


1 


Principal components of dispersion 


2 


Principal components of correlation 


3 


Principal factor analysis 


h 


Image analysis 


5 


Canonical factor analysis 


6 


Maximum likelihood factor analysis 


7 


Alpha factor analysis 


8 


Minres factor analysis 



Table III 

Factor Canonical Variate Rotation Codes 



Codes ' Type of Rotation 

1 Quartimax 

2 , Varimax 

3 ^ Equamax 
k Promax 

5 Multiple groups 

6 Orthogonal centroids 

7 Orthogonal bounds 

8 Orthogonal Procrustes 
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Table IV 

Power of .01 for Convergence Criterion 



, Power X 10 Criterion 

11 .0T9'*329 

12 .0630957 

13 .0501188 
lU ' . .0398107 

15 .0316227 

16 .0251189 

17 .0199526 

18 .01581*90 

19 .0125893 

20 .0100000 



Note: It so happens that the criterion is modular in the 
secon(L digit of the power of 0.1 while the first digit 
determnes the number of zeros after the decimal; i.e., 
use of 2h gives a criterion of .00398107 and gives 
a criterion of .0000398107. 
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Control Cards in Detail 



TITLE Card 



Cols. 

l-h 

5-80 



TITL these letters exactly 

tAny description that can be keypunched 

As mar^ TITLE cards may be used as the user 
sees fit. All cards must have the letters 
TITL in columns 1 to ^ 




NCONT 



Explanation 
PROB tnese letters exactly 

Number of variable format cards used 
(10 or fewer) 

The variable format cards are used to describe 
the scores as they appear on the data 
^ie. This information tells how many 
variable format cards are necessary to 
describe the format of thp scores. 

Number of contrast card ^sets 
(at least 2; at most 7) 

Each contrast card set describes a factor 
of an analysis of variance design except 
that one additional set is used to describe 
the continuous variables of the linear model. 
Thus a one-way analysis of variance model 
must have two contrast card sets: one for 
the design and one for the continuous vari- 
ables. A simple correlation problem must 
have a dummy contrast card set even if it is 
only one sample, and one set for tne con- 
tinuous variables. 

In many problems the contrast card sets will 
have only one card per set . 
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Cols 



10 



Codes 



INF028 



12 

17-18 



MPHINT 



KPRINT 



INFILE 



20 



IHF021 



INF09 



Explanation 

Number of individual sihnificance test 
control v?ards (at most C) 

For any tests calculated it is possible 

to enter ah individual significance tesc 
control card to cor^rol use of any of 
several expository techniques such as 
rotation, bat^tery reduction or the like. 
At most 6 tests may use these tech- 
niques, so that at most 6 individual 
significance test control cards may be 
used. 

Punching a 1 (one) forces-Tprinting of 
the cell means and variance; leaving 
a blank avoids calculation and printing 

Punching a 1 (one) forces printing of the 
reduced model matrix 

If the data _ire ou cards and submitted 
with the set-up cards, leave blank. 
. If the data are on a com^^er kept file, 
punch in the file number 

Punch a 1 (one) when only estimates of 
effects or regression coefficients are 
desired for all the tests. Otherwise 
leave blank 

This column controls what output is to 
be printed. Leaving" it blank prints 
the bare essentials <;r the test. 
Punching a } (one) or f (two) prints 
additional output, '^ee Table I 
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INFO10 
inF0l6 



This column also controls vhat output is 
to be printed. Gef* Table I 

This column controls ^egression in analysis 
of covariance and m' st^ be used for all 
covariate adjustments of data. Punching 
et 1 (one) produce*: the classical analysis 
of covariance adjustment. Punching a 2 
(two) produces a reduction of the hypothesis 
sum 'of squares by the error regression 
weights. Failure to enter a code results 
in the class i'^ai adjustm^^nts. 
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Cols Codes 



28 



INE035 



Ex-planation 

Punching a 1 in this column prevents copying 
of the data file for reanalysis. This is 
particularly useful for large data sets 
when no reanalysis is to "be done; it saves 
time on the machine. 
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3. Contrast Card f>etr> 

The description beL^w showr> ^hni il ^} rooalr^M inf- nv^^ . r '^'^ 
rriore than one card, hence th^^ ter-. '\"ir i .^'<'t" r'tt h^r *ri'ti: ""in. 

One set is required for the <^ontinuous viri'ibl^r 'ini -iv 
factor of an ANOVA desle;n. If the dnfi comprise 'i slnr.le ivu-ri* 'rz ofu-n 
happens in factor analysis or corre.ition analysis, 1'^ Is s'll^ 

^ ■Mil 

to include a contrast card set tor "factor of a iesirn ; ^n-* 
the dat m as if it were a one-wa^ lesip.n vlt:. a sin^;ie le/eJ. " -r *hir> it 
may be convenient to have a 1 at sur.e position in eac/ri ^lata r-c^r:; ' r 
any constant digit can be used and recoded, even a blanV. Ilie -vl-r of 
the card sets must be the same r;s the order in which \he d-ta f'l.- >r i^. 
appear on the variable format statement. 

It is not necessary that "far-tors" of a .design appear rl- 
significance test card. This makes it possible to u^e a "factor" index 

for file elating or in later reanaiyses of the data. Use of factor 
code f^r^lting should be limited sin^e the ^ulled cases are r-;ected 



and listed on the output. 
Cols Name 
1 NT/VBLE 



ICOiNT 



Explanation 

For a factor of an AIJOVA design: Any 
letter except V or V, but no two 
factors may use the same letter 

For the continuous variables: 
'Yhe letter V only 

Contrast codes for design variables 
= b Regular contrasts: Also called 
nominal and deviation contrasts 
Example: for t.hree levels of a factor 

mean 1/3 1/3 1/3 

df^ 2/3 -1/3 -1/3 

dfj -1/3 CV3 -1/^ 

= 1 Special contrasts. Any set of contrasts 
the user wants. Tlie matrix must be fed in 
as described below and must follow this card 
set immediately 
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Cols Name / Explanation 

= 2 Orthogonal polynomials (integral values 
only). This option generates orthogonal 
polynomial contrasts and design parameters as 
\ integral valuoo. The number of levels is 

limited to 19 for numerical reasons. The 
metric of the polynomial (i.e., the lir^^^r 
coefficients) must be supplied and must follow 
this card set immediately. 

= 3 Reverse Helmert contrasts: also 
called difference contrasts. These 
contrasts are as follows: 

mean 1/^ l/^ 
df, -110 0 
dfi -1/2 -1/2 1 0 
■ df^ -1/3 -1/3 -1/3 1 
It should be noted that the fractions 
involved in these contrasts are not 
exactly representable on computers. 
This may cause serious numerical problems 
due to rounding errors when the number of 
data cases is large, say a c6uple of 
hundred 

= 1( Special design parameters. Any set of 
design parameters the user wants. The 
'matrix' must be fed in as described below 
and must follow this card set immediately. 

= 5 Orthogonal polynomials (decimal values 
only). This option generates orthogonal 
polynomials with normalized coefficients. 
Coefficients can be generated up to ord-r 39, 
and the polynomlais of order 39 are accurate 
at least up to dagrec % The accuracy of the 
higher order polynomials is not guaranteed. 
The metric (i.e., the linear coefficients) 
must be supplied anvi musl follow this card 
!:,ot immediately 

For power functions of continuous variables 
= 0 no powers generated 

= N (after V). When the contrast card 
set describes the continuous variables, this 
number describes the number of powers of 
each variable that wi]l be generated to 
form a new variable; that is, if n = 3, 
the variables V, V? and v3 will be available 
for analysis. It wilJ also be assumed tha.t 
there are N times the number of variables 
(i.e., 3 X V) in the new analysis and that 
they are in the order V, V^, V^. 
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Cols Name 



Explanation 



U6 NLBV 



KRDNA 



Number of* levels for the factor or, after V, 
the number of variables (not includinr^ 
powers ) 

= b if the dummy or, after V, continuous 
variables are to be labeled by n^omber 

= 1 if the dummy or after V, continuous 
variables are to be labeled by names 
which will be Supplied immediately after 
this card set 
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KLINT 



11,12 
and 2 
column 
fields 



LEVSUB 



the next 
2 column 
field 



in 2 

column 

fields 



RECODE 
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For making linear transformations of the con- 
tinuous variables of an analysis. Three 
options are available. 

= 1 ordinary linear transformation. Coeffi- 
cients are expected in the original order of 
variables. [ / 

= 2 Potthoff and Roy transformation with the 
weighting matrix. The weighting matrix 
precedes the matrix of the transformations. 

= 3 Potthoff and Roy transformation with the 
inverse of the weighting matrix. The inverse 
of the weighting matrix precedes the matrix 
of the transformations. 

When the dummy or the continuous variables 
are to be subdivided into partitions, 
'this is n\Mber of variables in ejach 
partition. The number of degrees of 
freedom or variables are punched in 2 
column fields and must account for 
all degrees of freedom or the total 
nuiiber of variables. If no partitioning 
is done, ignore 

,, (two commas). This is necessary whether 
or not the variables are partitioned. If 
the variables are not partitioned, the 
commas go into columns 11 and 12 

Recoding of factor level identification or 
reordering of the variables 

For recoding the level identification, 
enter the code as it appears in the data 
file in the order in which the codes are 
to be remembered. For example, if the 
codes are 15, 5 and 31 and these are 
to be recoded as 3, 1 and 2 make the 
entries in columns 11 through l6 as b 53115 • 
Use two column fields throughout and account 
for all levels. Zero is an admissible level 
code. 

.32 
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Ccls 



Name 



Explanation 



in the 
next 2 
column 
field 



the next 
20 or 
less 
columns 



AHAME 



For reordering continuous variables, 
enter the serial numbers of the variables, 
as they appear on the data file, in the 
new order. Use two column fields through- 
out and account for the total number of 
variables. When several operations are 
done to the variables in sequence, the 
sequence of operations is 

1, raising to powers 

2, making linear transformations 

3, reordering 

1^ , partitioning 
5, naming 

In employing these features, the user must keep 
in mind such problems as (l) linear trans- 
formations must include all powers of 
the variables, (2) partitioning takes 
place on transformed and reordered variables 
etc. 

When variables are raised to powers, columns 
U to 6 must contain the number of original 
variables, but all subsequent operations 
must take account of the original variables 
and their powers. That is, if there are 5 
original variables and these are raised to 
the jWwsTd power, the program will expect 
15 Linear transformations, 15 variables to 
reoMer, Y5 variables to partition and 
15 names 

(two commas). This in necessary whether or 
not the varlabJcs are partitioned and 
whether or not r(^ordering or recoding has 
been used, IT neither feature is used, 
the commuG go into columns 13 and 1^ 

Any alphameric description of the factor or 
variable set. This may be as long as 20 
characters 
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Cols Name ? Exp3 anation 

Last . (a perlo'i) 

Coltimn 



h . Contrast Card Set Options 

^ese descxiptions govern use of options found on the contraat cards. 
Of the k options, use of one precludes use of the others except for the 
variable names option. When variable names are supplied, the naiues <"ar''is 
must precede cards used for other options; All options requirimr cards 
must follow the contrast card set to which they apply and precede any 
successive contrast card set, 

a. Variable names cards 



The use of this option enables printing of names on the output. 
Both continuous variables and design parameters (or contrasts) 
may be named, although dxmmy parameters for ANOVA interactions 
will always be numbered. 

Names 'will have 8 characters (including blanks / and 10 names 
may be put on a card. When dummy parameters are named there 
must be as many names as there are degrees of freedom. V/hen 
continuous variables are named there must be as many names 
as there are variables. When the continuous variables are 
reordered, the names must be supplied in the new order as the 
names are not reordered. When linear trans fomations are made, 
the names will be affixed to the transformed variables 
(including analysis of covariance). When powers of the variables 
are generated, the program expects as many names as there are 
variables times powers. 

b. Special contrast or special design cards 

This option enables use of special one-way contrasts and design 
parameters. They are entered as square matrices with as many 
rows and columns as there are levels of the factor. The first 
ccntrast or design parameter entered is always that for the 
constant term or grand mean and is the first row of the matrix. 
The weights for obtaining the coustant term or grand mean nee' I 
not be all equal nor do the contrast coefficients or design 
parameters need to sum to zero, although this is usually 
desirable. 
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/ 



The elements of the matrix are punched in 10 column fields, 
8 elements per card. Each row must begin on a new card and 
may continue on as many cards as necessary. 

c. Orthogonal polynomial metric 

The coefficients of the linear polynomial are called the metric 
and are usually (depending on the expariment) the integers 
from 1 to N indicating equal spacing of the levels of treatment. 
Equal spacing is not necessary, however, and the metric need 
not be successive integers. 

The coefficients of the linear polynomial are ente/ed in 10 
column fields, 8 per card and continued on as many cards as 
necessary. » 

d. Linear transformation of the continuous variables and Potthoff and Roy models. 

This option allows for testing linear combinations of the 
continuous variables. To use option 1 it is necessary to use 
all the variables in the regression equations and to have as many 
regression eauations as there are variables even though some of 
the coefficients might be zero or one. The regression equations 
are entered as a transformation matrix, with the equations as 
rows. Each row begins on a new card, the coeffic^ ts are punched 
in 10 column fields, 3 per card, and on as many c as as necessary. 

When using option 2 or 3 for Potthoff and Roy models, the weighting 
matrix or its inverse precedes the matrix of the transformations 
The elements weighting matrix or its inverse are entered a row at 
a time, punched in 10 column fields, 8 per card the same as the_ 
transformation matrix. It is essential that the weighting matrix 
precede the transformation matrix. 

'^en reordering is used in conjunction with linear transformations, 
the reordering takes place on the transformed scores. 

5'. Significance Test Card 

The significance test card describes the model under which the data 
are to oe analyzed ani tested. Linear models in three styles of calculation 
are available: analysis of variance, correlation, and factor analysis. 
The calculations be made are indicated by the acronyms MOVA, CORREL, 
or FACTOR in the lirst few cola-„ns of the significance test card. 
Calculating interaction tables in a factorial sampling design can be 
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indicated by the acronym IIITABL. Variance reduction studies can be 

calculated by using the acronym VARED. 

Use of the acronym ATiOVA invokes solution of the linear model as an 

analysis of variance problem. Either one variable or many may be analyzed 

in the model. The null hypothesis test utiMzes Wilks* iaxnbda 

criterion as approximated by the F distributi >n. Although the lambda 

distribution is ostensibly multivariate, its degenerate cases, uni- 

2 

variate ANOVA, Hotelling's T , Mahalanobis* distance D, discriminant 
analysis and "Student's" t are aJ.1 handled automatically. 

Use of the acronym CORREL invokes solution of the model as a corre- 
lation problem. The null hypothesis test is Wilks' lambda criterion as 
is the case for AlfOVA problems. Again, the degenerate cases of canonical 
correlation: multiple correlation, product moment correlation, biserial 
correlation and point biserial correlation are all handled automatically. 

Use of the acronyp SERCOR invokes a serial correlation analysis 
model. Serial correlation is used to analyze characteristics of 
observations within samples (in contrast to ANOVA which analyzes 
differences between samples). Serial correlation models are a mixture 
of ANOVA and CORREL models. To calculate a biserial correlation between 
a dichotomous characteristic. A, and variables, V, the model statement 
SERCOR:W=0;A=V. 

is appropriate. This states that the effect for the grand mean, W, 
is swept out of A and V because W is to the left of the semicolon. 
The effect A is not swept out of the variables, V, «o that A = V is 
analyzed as a correlation model. This is a simple model which could 
also be calculated as either an ANOVA or CORREL model. 
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More complex serial correlation models are coinmon for a model with 
several samples. A, and one characteristic, B, common over the samples 
the following model is appropriate 

SERCOR.-W=0,A=0;B=V,AB=V . ^ 
Here W and A are to the left of the semicolon to remove the effects of 
grand mean and between sample variation from the data before analyzing 
the characteristic B. The effects for B are not removed from AB and V 
because they are to the right of the semicolon; likewise AB is not removed 
from V. The effect AB=V is a test to check the homogeneity of regression 
of B onto V in the various samples of A, Further discussion of the uses 
of serial correlation is found in the reference cited in the bibliography. 
Serial correlations can be performed with all the variations of the linear 
model: covariates, rotation of canonical variartes, Potthoff and Roy models, 
dimension reduction, etc 

Use of the aconym FACTOR invokes a factor-analytic decomposition of 
th^ data. Gevex-al types of solutions are available as displayed in 
Table II. The acronym FACTOR must be immediately followed by a digit 
froin Table II to (denote which factor decomposition is to be used. It 
is possible to denote a rotation procedure by following this digit with 
a number sign i fr) and another digit to denote type of rotation, vmen 
cominunaiity procedures other than squared multiple correlations are used, 
the individual significance test card also must be used. When less 
^hari all the factors are to be extracted, the individual significance 
test card must be used. 
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Use of the acronym INTABL generates interaction tables foi^ ANOVA 
problems. Interaction tables are generated as a nesting procedure with- 
out including the grand mean in the model. Thus the statement ,WAI3=V1, 
generates all the means of the variables in V for the AB interaction 
of an AHOVA model. It is possible to generate many interaction tables 
in the same model statement. The estimates used for constructing ^ 
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interaction tables are the raw estimates and are not subject to analysis 
of ccvariance adjustments. Interaction tables adjusted for covariates y 
may be obtained by u: ing regression transformations on the original 
variables. 

Use of the acronym VARED generates a study of the sum of squares 
for hypothesis before and after reduction of the mode^ to orthogonality. 
Mode, statements follow the rules for AI.'OVA and CORREL models. VARED 
studies are not subject to analysis of covariance adjustments except 
by using regression transformations on the original variables. 

On the remainder of the card, special symbols, letters and numbers 
are usea to designate the T.odel to be analyzed. 

Letters . The letters used are those from the contrast card set. 
Iliat is, use oftne^etter A assumes that there is an effect and a factor 
jf the design to be called A and that there is a contrast card set that 
his "A" pur.che.d in column 1. 

".Tiere is also a contrast card set for the continuous variables which 
has a "V" in coliinn 1. '^len the phrase ,A=V, is punched into the card, 

ir of 3a,.-trvr. for t^e A effect will bo generated and tested for its j 

/ 

r-f-rescion or. ti.e v-i-'l 'ill es V. 

If the desi^:n U a factor' ai, there nay aiso be a B effect to test. 
r,r this we punch the phrase ,3=V, ir. the significance test card. To 
t"3' the interaction c ffect we^ punch the phrase ,AB=V, in the card, 
rh- une of two or more letters adjacent generates the Kronecker product 
,,r -he main e'Tect .iesir;n parameters necessaiy for testing interactions. 

^rder of Kronecker product . 'iTiis prograjn generates Kronecker products 
in thf' order letermined by the order of the contrast card sets and 
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ultimately by the input data record. Consider the term Ah, I:" t.^c 

contrast cards ^et for f-ictor A precedes that for f'lctor (whi^. ii: 

say the identification code for factor A is to the left of M^.o i h-nti- 

ficEtion code for factor P on the data file) the AB durnmy parameter.-- are 

generated and the estimates printed in the order ^^^i^ ^J^l' ^-^^i^***' 

a b, , a,b^, a^b^,...,a b^^a^b^,, ,,,a b^ . The same order is gc^nerated 
nl' 12' 22' ^n2'15' 'nm 

whetAer the interaction term is written BA or AF. 

Equal sigrx . An equal sign ( = ) is used to separate the hypotliesis 
variables from the error variables. ITie hypothesis variables are 
designated by letters and numbers to the left of the equal sign; the 
-error variables are designated by letters and numbers to the rip;ht of 
the equal sign. \^ 

Commas, colons, and periods . Commas (,) are used to separate the 
tests from each other. A colon ( : ) is used to separate the model 
acronym frpm the te<=ts; and a period (J is always the last character 
in a model statement. 

Grand mean and zero . The constant term or grand mean may be included 
in the model by using the phrase ,W=0, where W indicates the constant 
term as a hypothesis- and 0 (a zero) indicates that there are no error 
variables for the hypothesis and no test to be made. The use of zero 
as an error terra is the way in which a set of variables may be included 
in the model as if they were hypothesis variables but not tested. 
^Jt\en 0 is used as a hypothesis, it indicates that the error term is 
to be subjected to factor analysis. 

With this information it ib possible to write a significance 
test card for a simple factorial analysis of variance: 
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(a) ANOVA:W=0,A=V,B=V,AB=\ . 

'I sif";ni f icmee test cin'l for x C'ictr^r 'UKtlyr^is: 

(b) FACTOHl:W-y^,^ V. 

::;anibered partitions , l^e contrast cirJ optica for partitions mak s 
it po^'sible to subdivide tne f.els Ito.^^r. par-imeters and continuous 
variables into several subsets. Tbt- iMrtitions are written on the 
si.^.nificance test card as Al , A.\ or \I, V3, hi, etc., the number 
referrin?:; naturally to the oidinal position of the partition (this use 
jf niimbers can be easily ^vnifusei vLtr. the use of nu.T.bers in n^^sted 
rr.oaels ) . 

With this inforr.ation it ir. poor.iblt^ to write the following models 

'ir. 1 raany others . 

(c) i^artiticn^M anal%-i3 of variance (such as orthogonal 
po 1 y non 1 al e t ) 

A:1OVA:W=0,AI=V,A =V,B=V,AIB=V,A2B=V. 

^ i) i'^^rrel.at i r.oda^^ 

, , :!or.' nel ty r. in a one-way design 

)VA : W=^ , A^y^ , V 1 '-'^ , A7 . = VP . 
(f) A princip<a. faol-r .uiaiysis or. several samples 

FkC?'^^ >.\: /.A=/./ V. 
ir) Correlati( ri analyias from several samples 
:'OBHrL:W=!2i,A=0,Vi=V. * 
, ^[^^) ^iri. l> '^i>*^i t" ^> several effects to produce 

: i- .'^i'xres. In an anaiysi" of varKince with nev-ral factors, one 
: '.,.-r .'-i^^r int^-r-Lcti wn.- and ♦'^r.t then ir.ultanenu-.ly ''^^r Uie 
r, - I v't>';ab:t*s Vi+V^ 
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Fur instance, 

(h) ,A!^+A^+^^+AP^=v^ , 

When orthoi^onai poIynomi'iL' :in> part It i')r/^i, '^hir -i: ro -tl lov*^ - 
r^ooiing of high oni'-^ polynomiaJ for oimaltrinoour tcr.tinr, 'xc ! - 1 ! v:-., 

(i) ,Pi=v,p:>=v,P3+r^i-^r^;=v. 

\^en + is used belweer' two sets of variables the variables are 
pooled before any other operation takes place: For ex.ample VI /v: +V 
indicates that both V? and V3 ar. mtly oovariates for VI. 

Wlien the parameters for a give, effect are used both singly 
in some tests and pooled in other tests, the order in which the 
parameters are presented on the significance test card must alwayr 
be the 'same for every test in which the parameters are used. 

W and numbers for nested analyses . The letter W is commonly used 

in statistical literature to indicate nesting; BWA indicates tha* 

seversLl samples, B, are nested in each of the levels or samples of A. 

This notation is expanded slightly here as follows. 

(j) ,BWA1=V, indicates a test of B within the first 
level of A. 

(k) ,BWA1+PWA2=V, indicates a pooled test of B within 
only the first two levels of A, 

(i) ,BWA=V, indicates a test of B pooled for all levels 
of A, 

(m) ,BWAC, indicates a test -f B pooled over all cells ^ 
of a two-factor design. A"". 

(n) ,V1WA1=V2, indicates a to^t cf the correlation between 
VI and V2 within the fir^t level of factor A. 

Here, a number, used after a letter which follows a W, indiratoG 

the level, not a partition. A number used after a letter but before 

a W (or in the absence of a W) indicates a partition. 

4'^ 
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Slash (/) for analysis of covariance . The use of a slash indicates 

that all variates after the slash and before the next comma (or colon 

or period or minus) are to be covariates for the test indicated. A 

slash and variables immediately following the acronym indicates that 

all the tests in ttie model have ♦.hose variables as covariates. This 

makes possible tht fo] lowing types of models. 

V.) Analysis of covariance in a factorial analysis of 
variance 

ANOVA:W=0,A=V1/V2,B=V1/V2,AB=V1/V2. or 
AiI0VA/V2 :W=^,A=V1 ,B=V1 ,AB=V1 . 
(p) Partial correlation 

C0RREL:W=^^,V1=V2/V3. 
(q) Principal components analysis with covariates 

FACTOR1/V2:W=0,0=V1. 
(r) AIJOVA:W=0,A=V1,B=V1/V2,AB=V1/V2+V3. 
Ampersand (&) for extension . The use of an ampersand indicates 
th^t ail the variables after the ampersand and before the next comma 
' r - .on - leriod cr minus) ire to be used as extension variables 
t: nne r^anoni -'il variates of the test indicated. An ampersand immediately 
foiiowin^, the 'icronym inaicates that extension is to be done to all 
th^ canonical variates in <xLl the tests in the model statement. Thii> 
"takes :t .)>^ssible to write the following models. 

vs) lasrriininant analysis of VI and A extended to V2 
ArjOVA:W=(^,A=Vl&V?. or 
/aJOVA&V:^:W=0,A=Vl . 
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(t) C.'tnoni.'u rrelativn l-'Wf/en F und VI ^-xt^rrK-i to 

^\)HHEL:V=^^,B=Vi&V. + 
(u) Faotor 'irrilyii' - Vi t;xtondo<i lo VP. 
^ FACTHfJi :W=0,0=VI^V , 
I-jumber sign f/ff) for r.Mdtion . 7h*> use of after the a^-ronyn 
FACTOR and its iir,it iniicites *ht* faetors are lo be rotate. i. A 

digit must follow to ienoto which rotation scheme is to be isod. 

/ 

Table III lists the typer -f r< tations available, 'lliis mfikes it possible 
to writ*^ the model. 

(v) Alpha I'aetor rinaiysir» with equamax rotation 
FACTOH7/M:Wr:^,0=V. 

(Fetation of canonical variates in AJ^OVA or COPREL 
can be done uniy by usinp. the individual sir/iif icanc^ 
test card. ) 

Asterisk (^) for coT np< ments cj f variance . In many components of 
variance models, the error term for a test can be generated as a 
regression siim ^f squares. More simply, the sujii of squares for hypothesis 
of one test may be the sun of squares for errors of another test. 
Therefoi-e, to make it easy to cai -dilate the sum of squares for error 
we use the notation A*V denote ihe regression of cont-^nuous variables 

V on dummy paraineters A tn fom an error term. This notation is allowed 
only on the right of an equal si^^n where errors are designated. 

It should be noted that this yru^rrun will not solve multivariate 
co^inonents of variance models wherr* the numb^^r of continuous variables 

V is more than the niimber of lummy parameters A or degrees of freedom 
f r- ro r , 
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Minus (-) for dimension reduction . In multivariate models, it may 
be desirable to remove the significant canonical variates of one effect 
from the variables to be analyzed for another effect: to wit, the 
significant canonical variates of an MOVA interaction rc lOved from con- 
sideration in a test of a main effect. If A and B are main effects in 
an ANOVA design and AB their interaction the statement 
,A=V-AB=V . 

will generate a test for A independent of any significant interaction 
effects • 

Semicolon ( ; ) for limiting sveeps in serial correlation. In ANOVA 
and CORREL models all hypothesis terms in the model statement are swept 
out sequentially from left to right leaving a residual sum-of-squares of 
the variables} In serial correlation, the hypothesis terms which lie to 
the left of the semicolon are swept out but those hypothesis terms to the 
right of the semicolon are not swept. The assumption being that liypothes 
statements to the right of the semicolon refer to characteristics of the 
observations in the samples or populations. 
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Special notes * 

(a) When effects are partitioned and repooled (for example V1+V3) 
the pooling must always be presented on the significance test card 

in ^he same order. That is, you cannot state A=V1+V3 and B=V3+V1 
in the same model- You must state A=V1+V3 and B=V1+V3. Also, when 
effects are partitioned, they must stay partitioned: i.e., A=V is not 
an alternative for A1+A2=V when A has been partitioned. 

(b) In one model statement, a set of parameters may be used only 
for a hypothesis, an error, an extension or a covariate but not for 
two of these. 

6. Variable R)rmat Cards 

These cards describe the way in which the data on the observations 

appear on the data file. There may be as many as 10 cards to describe 

the format. The use of the variable format follows the customary 

FORTRA!J restrictions. 

Suppose the data appeared as follows: 

Cols. 1-5 information to be ignored 

6 level number :^f the first factor 

7 information to be ignored 

level identification of the second factor 
2^0^13 a datum on the first continuous variable 

3^l^_.l6 a datum on the second continuous variable 

19^21 information to be ignored 

22^25 a datum on the third continuous variable 

26- to end information to be ignored 



■y 

The information on these records can be read with the following 

variable foi-mat statement 

(5X,Il,lX,l2,Fl*.o;F5.0,3X,Flt.O) 
It will be noticed that all identification numbers are expected in 
fixed (I, integer) fomat while all continuous variable scores are 
expected in floating (P, decimal) format. 

About factor identification . The reading of factor level identification 
must precede the reading of variable scores. In the above e.cample the 
data record is arranged so that this occurs naturally. It is possible 
that the factor level identification is interspersed among the variable 
s:ores; in this case the "tab" feature of format statements may be used to 
read the factor level identification before the variable scores. If this 
cannot be done the data will r.eed to be rearranged to put level codes first 
in the records. 

All records must contain factor level information. When single 
samples are analyzed, it will be necessary either to include a constant 
on the file or to fake a factor with one level. This may be done by 
reading a number off the record and recoding it to 1 or b^ reading a 
blank on the record and recoding it to 1 . ^ 

On :BM machines it will also be necessary to note the following 
comment about reading integer variable scores in F format: for example 
reading the score 12 as F2.0. 

The score 12 read as F2.0 occupies 2 characters. When the program 
copies this score to save it for reanalysis, it copies 12.0 in two 
characters 2.0 which does not include the 1: The copy is overflowed. 
When the copy file is read for reanalysis the program will register an 
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overflow in subroutine DATVEC. To prevent this disability either 
manufacture the original file as 012 and read it as Fi.O or manufacture 
the original file as 12. and read^ it as F3.0. 
7. Data 



on 



The data may be on any file as long as the file is designated 
the problem card by a two-digit number in columns 17 and 18. Blanks 
in these columns indicate the data are on curds and follow the variable 
f o rma t s t at eme n t • 

If the data are on cards, a blank data record must follow the data: 

there must be as many blank cards as there are cards in the data record 

/ 

of ofie observation. If the data are on a tape or disk file an ordinary 
end-|of-file mark will do. 

''as a rule it is simpler to arrange factor level identification codes 
first on the data record before the continuous variable scores because 
these codes must be read in first in the record. This expedient is 
not a necessity on many machines because of the *tab " feature of 
FORTRAN IV compilers. 

8. Individual Significance ^est Control Card 

The individual significance test cards are used to control procedures 
which can be applied to particular significance tests and not to others. 
Iliese procedures are expository in nature and not usually subject to 
statistical testing. Arbitrarily, the number of statistical tests 

in the model to which these procedures can be applied is limited to six. 
Each individual significance test control card is presented with 

the optionii which apply to it. The control cards are presented in the 

reverse order from the order of the tests on the significance test 

card (i.e., last to first). 
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Column Name 
1-2 ^'OW 



3-1. 



CRIT 



T-8 



MJVXFAC 



10 



Information 

The number of the significance test to 
which this information applies. To 
obtain this number from the significance 
test card start at the period and, pro- 
ceeding backwards toward the acronym, count 
the number of equal signs u.itil the test 
which applies is reached, use that number 
in columns 1 and 2. ^ 

The probability, multiplied by 100, used to 
control (1) additions and deletions of 
variables in battery reduction, (2) 
dimension reduction (use of minus signs) 
in ANOVA and CORREL and (3) the number of 
canonical variates to rotate in ANOVA 
or CORREL. 

Type of battery reduction procedure for error 
variables 

=1 Efroymson's stepwise addition/elimination 
procedure 

=2 Wherry-Doolittle stepwise addition 
procedure 

=3 Step-up variable deletion procedure 

Number of factors to be utilized in factor 
analysis tests, rotations and extensions. 
This does not apply to "the two principal 
components solutions as they always obtain 
as many components as there are variables. 
It does not apply to principal factor 
solutions as it always obtains one less factor 
than there are variables. All other procedures 
require an entry. 

Statistical test procedure for factor analysis 
=0 none 

=1 Rao's test, for ct-ronical factor analysis 
only 

=2 Rippe's test 

Type of rotation technique for factor analysis 
or multivariate ANOVA or CORREL. (Use of 
this is not necessary if the H (number sign) 
is used after the acronym FACTOR. ) See Table 
III for the usable codes. 
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Column Name Information 

ll^ For direct Or indirect rotations 

=1 Direct dotation of canonical variates 
or factor analysis rotation 
=2 Indirect rotation of canonical variates 

For taxonomy of ro\ation 

=1 for taxonomy ofv variables 

=2 for taxonomy of Vroups (in MAIIOVA only) 

I'j^lQ Criterion for uniquenessV covergence in canonical 

factor analysis. This\ criterion is entered as 
an exponent (multiplied\by 10) of 0.1 so as ^ 
to give a wide range of Values to convergence. 
The entry may be any two-^igit number (see 
Table IV). For example, the entry 31 
produces a convergence criterion of 0.1 * = 
.00079^329. 
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Communalities supplied for factor analysis 
=1 if squared communalities are supplied 
for factor anal^/sis. Only three procedures 
use this option: principal factor analysis, 
canonical factor analysis and image analysis. 
If communalities are supplied, they must 
immediately follow this card. 

=0 if squared communalities are to be the 
squared multiple correlations between the 
variable and all other variables. IT IS NOT 
> NECESSARY to use an individual significance 

test control card to have squared multiple 
correlations used as communalities. All 
common vari.:^.ce factor analysis procedures 
have an automatic default to squared 
multiple correlation. 

=1 if weights are supplied for multiple groups, 
orthogonal centroids or orthogonal bounds 
rotations. If weip.hts are supplied, thoy 
must immediately follow this card or the 
supplied squared communalities, 

9 . Individual Significance Test Contr o l Card Options 
a. Supplied communalities 

Communalities (squared) must be entered 8 per card, in 10 
column fields, and in the order in which the variables are to be 
analyzed. 

o ' so 
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b. Weights for rotations 

The weights on each variable which determine a transformation 
vector are entered 8 to a card on as many ca.ds as necessary. Each 
weight is entered in a 10 colunrn field. In each transformation 
vector there must be as many weights as there are variables and 
each transformation vector must begin on a new card. 

The number of transformation vectors to be entered is determined 
by the problem. For AWOVA and CORREL problems, there must be as. 
many vectors as there are canonical variates; i.e., the minimum 
of the number of variables (degrees of freedom) for the iiypothesis 
or the number of error variables. Fof factor analysis the number of 
sets of weights must be either the number of variables or the number 
of factors specified in columns T and 8 of the individual significance 
te?)t control card. 
10. TITLE Cards (optional) 

TITLE cards may be used to identify reanalyses. 
!(.. Reanalysis Card (optional) 

This card allows, in the same run, reanalysis" of the data. Two major 
changes are allowed in reanalysis: (a) change of kny or all contrast card 
sets and (b) change of the model statement. There are some features of the 
first .analysis which cannot be changed on reanalysis: (a) the number of 
contrast card se^s used cannot be increased (although some factors of the 
design may be ignored in the mode: statement); (b) the number of levels 
or variables may not be increased (although some levels or variables 
may be ignored on th^ model statement); and (c) the format of the data 
(the variable format card) may not be changed. 
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Cols Information 

ANLY, these letters exactly 

6 If a new significance test card is used to 

alter the model analyzed, punch a 1 (one) 
in column 6. (Each change of model must 
be accompanied by a reanalysis card) 

Q The number of changed contrast card sets 

which are to be used in reanalysis. This 
feature can be used to alter several features 
of the model. 

(a) The contrast used for a given 
"factor" of the design may be changed. 
The letter designation of a "factor" 
may not be changed. 

(b) The way in which the dummy parameters 
and continuous variables are parti- 
tioned into subsets may be changed. 

(c) The names of variables may be changed. 

(d) The transformations of variables may 
be included or changed. 

To change these features simply follow 
the original instructions above. 

The following- features are the same as described for the PROB card 

and appear in the same card columns. 

10 

^2 Print cell means and variances (1 means yes) 

Print reduced model matrix (1 means yes) 

Print only estimates or raw regression 
coefficients (l means yes) 

22 Controls printed output. See Table I 

2^ Controls printed output. See Table I 

Type of covariance adjustment (O or 1 does 
^ the classical adjustinent, 2 uses the error 

regression weights for reduction) 



The number of individual significance test 
control cards for this reanalysis 
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The options of columns 10 to 26 must be reinstated for each reanalysis 
because doing a reanalysis erases the controls of the previous analysis. 

12. Several Problems 

Many sets of data may be analyzed in one run. The cards for each set 

may follow each other. 

13. FINISH Card -^^ 
A card with FINISH in columns 1-6 may follow the last problem. 



PrograjiuDing Nott.s 

The VARAN program is written in FORTRAN IV for an IBM 360-65 computer. 
Insofar as possible it is written to be compatible with any computer which 
use3 FORTHAIJ-like compilers. Such features as may nee^'^ changing are 
accessible without much labor. The following list shows some of these 
feature's . 

1. Ml multiple precision statements are titled "DOUBLE PRECISION." 

2. All input-output units are designated hy integer constants which 
can be changed in the main program. 

3. All dimensioning constants for changing the size of arrays are 
centralized in tne main program and four executive subroutines. 
This allows for easy changing of size parameters and adaption of 
the program to special problems. The arrays and their dimensioning 
constanLs are listed in Tables VI through VIII. The program is 
distributed in three sizes, Tiny, Standard, and Large. The FORTRAI] 
ieckr> are the standard size while instructions are on the tape 

for implementing the tiny and large sizes. 
The program is put together in two basic sections, one for the clerical 
part of setting up the model, and another for doing the mathematics of the 
ooiution. .The clerical section is controlled by the executive subroutines 
nUBr:XJ , 3U&EX3 and SUBEXli. The mathematical section is controlled by the 
e'::ecutive subroutine SU3EX2 with tnt. factor analysis subsection controlled 
by SUBEX5. Table V lists each of the subroutines with a brief description 

.r itr. T)urDose. This table lists each major executive subroutine and after 

i 

it, the r.iscellaneous subroutines used oy it as well as the major subroutine 



which are called by the executive subroutine. The order of che listing 
ro^jghiy dictated by the rollowing considerations: 

1. the order in which the program executes 

2. the order in which the overlay structure is put together, and 

3. the o> ' v;hich the FORTRAN decks are listed on the mailer 
tape (the -)\cepcion here is that th^ executive subroutines are 
at the heat of the mailer tape) 

Comment cards quite liberally to explain the functions of 

subroutines. A' ^-^^ iiead of each subroutine there are several comment cards 
which describe the function of the subroutine and how it operates. 

Accuracy of calculations > The programming sequence was arranged in such 
a way as to be convenient and accurate for multivariate problems. Therefore 
some of the univariate side statistics have been manipulated excessively 
outside the main flow of calculations and do noc serve as checks for 
nunerical accuracy. Accuracy checks should be made on the multivariate 
calculations . 
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A...er..blv lan,Tua.-e subroutines . ' Tvo of the subroutines MOVCHR and BINBCD 
are w?'itten in assenbler code because FORTRAII IV does not perform the necessary 
operations. Another IFRI-IA is word size dependent. Listings of these sub- 
routines are provided from the mailer tape but brief descriptions are also 
included here. 

MOVCHR. The calling statement is: Call MOVCHR (NS,SA,NR,RA ) . The 
routine moves the NSth character in 3A to the NEth position in SA. In the 
calling sequence NS is the integer position of the character to be moved, SA 
is the address of the source vector; UR is the integer position to receive 
the character, RA is the address to the receiving vector. 

EXAMPLE: 

Dimension hih) , B{li) 
where A contains A-BCDEFGHIJKLMNOP and B contains QRSTUVrfXYZ123!*56 

Call M0VCHR(12,A,3,B) 
Hesult: A unchanged 

i< contains QRLTUVWXYZ123!^56 
PIIIBCD. This is a real function which takes an integer, I, and converts 
it to alpha fonr.at, right .Justified. The following exfmple is given in IBM 
hexadecimal notation. 

Example X = BINBCD( 1 ,1'^) 

Input '^^^P'-'^ 
- ^ omnorni X = i^oUoljOFl 

i = noooooor x = i.oUoitOFfi 

Alth3u,-h the a ^hors are not c-rtain (having fallible memories) to the best 
.f ^ar knwwl.>d,:e VARAIJ uses only a single precision space filling option; 
tr.ir, is wny the secona argument is in in VARAII . The routine as originally 
ievi.ed is .-enrralized to allow g.^neral integer conversion, zero fill, space 
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fill, insertion of implied decimal point, and uses sin^^le or double precision 
results. 

IFRMA. This is a function which converts integral numbers in AIPHA 
format to integral numbers in INTEGER format. Zero is returned if a space 
is input 

Examples in HEX(36o) code: 



K = IFRMA{X) 



Input 



Output 



X = FlliOl^OUO 
X = F6U0U0I1O 
X = HoliOUoliO 



K = 00000001 
K = 00000006 
K = 00000000 



Table V 



i^ub routine Nanie 



Function 



Called from 



MAIN 



Main program 



IvlATA 



Aid matrix A to matrix B and put in C 



MANY 



M/\TX 

Miltiply matrix A tines matrix B and put in C 

'w>ites out error messages 
IMPIIIV 

Subroutine to obtain an improved inverse 



MAilY 



MANY 



MANY 



i:iv: 



Inversion of asymmetric matrix-Gauss pivotal 



IMPINV 
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IFRMA 

Function convert alphabetic characters to integers 

MOVCHR 

Moves a character 

LOOKUP 

L^earch table for value equal to X 
BC^BCD 

Converts binary to EBCDIC (IBM 360) 
JUBEXl 

Executive subroutine to set up data vectors 

RDCOiiT 

Head contrast cards and contrast card options and set up one way 
design parameter matrices 

DFVDES 

Obtains leviation contrast design parameters 

I'.PCON 

Read? in special contrasts (transposed) and computes design 
variables for special contrasts 

OPTPOL 

Computes orthogonal polynomial design variables 
HELDE" 

Obtains reverse Helmert design parameters 

58 



MANY 



mNY 



MANY 



RDSIGT 



mm 



RDCOIIT 



RDCONT 



RDCONT 



n">co:;T 



ibnvjtine Name 



Table V (Cont'd) 
Function 



SPDES 

Reads in special one-way design maLiices 
POTROY 

Reads in weighting matrices and transformations for 
Potthoff and Roy models 

HDGIGT 

Read the significance test card for the first time and set up the 
parameter codes for the terms listed* The scanning is done 
from left to right stopping at the period 

DAT^/EC 

Bet up the data vector for each data case 



Called from 
RDCONT 

RDCONT 
SUBEXl 



NEST 



(Generates an identity matrix for nested effects 
REMOD 

Print reduced model matrix- Generate data vectors for cells which 
are present 

aTTAT 

Finds the cell means and STD DEVS 

rUBEX3 

Executive subroutine to obtain cross products matrix and sweep it 



Accumulate and store cross-product matrices of the data vectors 



Perform Gaussian sweep by matrix blocks 

ZUhEXh 

Scan the significance test card -d set up for each individual test, 
Tirst scan up to the coign to get permanent covariates, extensions 
and rotation types which hold for all the tests* Then skip to the 
period and scan backwards- one test at a time» between commas and 
minus signs* 

HYERCO ^ ^ . 

Set up matrices for a hypothesis-error only, hypcthesis-cov, hypothesis- 
extension, or hypothesis-vector product analysis 

inHYCO 

Get up matrices for no hypothesis-error, no hypothesis-covar late, no 
hypothesis-extension, or no hypothesis /ector product analysis 



SUBEXl 
DATVEC 
SUBEXl 

SUBEXl 
MAIN 
SUBEX3 
SUBEX3 
MAIN 



SUBEXU 



SUB£Xi4 



INTER 



SUBEXl^ 



;et up matrices for interaction tab?es 
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Subroutine Name 



APLB 



Table V (Cont'd) 
Function 



Forin matrices for the analysis subroutines 

Reconstruction of matrices to conform to each other m terms of 
what has been swept 

GETM 

Get the I , J matrix from the SWP tape 

Subroutine for controlling analytic sequence 

ATOBT . , 

Subroutine to move matrix A to B and transpose it 



MAB 



Multiply matrix A times matrix B and put in C 
Multiply matrix A times matrix B-transposed and put in C 



M/\TB 



Multiply matrix A-transpose times matrix B and put in C 



OUTAB 

Subroutine to print out tables 



Called from 
MANY 

APLB 

APLB 
MAIN 



' Subroutine to get probabilities in several variance distributions 
including F, t, normal deviates and chi-square 

Subroutine to centralize IBM BRR3ET controls 

''''' Calculate the overall F-test for Wilks lainbda criterion, for _ 

multiple roots calculate chi-square for Bartlett's dimensionality 
reduction 



MANY 



MANY 



MANY 



MANY 



MANY 



MANY 



MJU^Y 



S>ubroutine to obtain eigenvectors and values of a matrix 



TRIDI 



MANY 



HOW 



rri-diagonalir.ation subroutim^ DWM I^^IT-UB 
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GO 



Table V (ContM) 
Function 



Called from 



hOW 



:et up simultaneous equations for eigenvector with eigenvalue E 



rlger.vilue i^ubroutine for tri-diagonal matrices DWM 151T-UB 



Subroutine :.o provi<ie estimates only or interaction tables 

. TMRED 

r-ubroutine to perform dimension reduction between separate hypotheses 
of MIOVA 

;OJHYP 

Subroutine to adjust hypothesis variables suni of squares and 
hypothesis-errors cross products 

^oJEPR 

Subroutine to adjust error sum of squares and raw estimates for 
covariates 

HEDHYP 

Subroutine to reduce hypothesis-errors cross products by error 
regression weights 

Er.TADJ 

Subroutine to obtain raw estimates adjusted for covariates from 
raw cross products 



ANAI 



Performs calculations for ATIOVA and CORREL 



DECOM? ^, 

Subroutine t6 decompose a symmetric matrix A into where T is 

triangular 

Subroutine to print AI^OVA calculations from ANAL 



PRCOR 



Subroutine to print correlation calculations from ANAL 

PRSER .14.. 

Subroutine to print serial correlation calculations from 

ANAL 

SUBEX^^ 

Subroutine to sort out factor analysis procedures 



TKAHV 



HOW 



SUBEX2 



SUBEX2 



SUBEX2 



SUBEX2 



SUBEX2 



SUBEX2 



SUBEX2 



ANAL 



SUBEX2 



SUBEX2 



SUBEX2 



SUBEX2 



MANY 



Subroutine to follow immediately after DETER 
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Table V (Cont'd) 

Subroutine Name 
DETER 

Obtain determinant of A 
JIMEQ 

Solve {A + PHI*I)X = B . A is symmetric • 

^^^'%ive- (A + PHI*I)V = X subject to V(TR)*V = CON (usually = 1.) 
Restricted least squares. 



:rfac 

3ubr 



'catine to generate first .five solutions of factor analysis 



Subroutine for alpha factor analysis: Michael Browne, circs 



L968 



MALFAl 



Maxinum likelihood factor analysis. Version 2. Variables giving 
Heywood cases are partialled out. 



MINRES 



Least squares factor analysis. Fixed unit weights. 
PRIFAC 

Subroutine to print factor analysis results 



RAO 



Subroutine to compute Rao's test of significance for canonical 
factor analysis 



RIPPE 



Subroutine to compute Rippe's test for completeness of factorization 



ROTCOn 

Subroutine for controlling rotations 



VARMAX 



utine tn perform quartimax, varimax, and equamax rotations 



Subro 



PPOiMAX 

Performs promax rotations 



WATCEN 



Subroutine for doing three rotations, multiple groups, orthogonal 
centroids and orthogonal bounds 



Called from 
MANY 

MANY 

MANY 

SUBEX5 
SUBEX5 



3UBEX5 



SUBEX5 



SUBEX5 



PRIFAC 



PRIFAC 



iJUBEX? 



SUBEX2 



SUBEX^ 



SUBEX2 
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'ub routine Name, 



Table V (Cont'd) 
Function 



F'actcr to lower triangle and orthogonal matrix 

PR I ROT 

Subroutine for printing results of rotations 

rXTEIID 

Subroutine to perform and test extensions 
HATRED 

Subroutine to do battery reduction techniques on error variables 



ELIM 



Subroutine to partial out variables from a canonical correlation 
problem 

Subroutine to calculate loss of variance due to orthogonalization 
of design or o'f covariance analysis 

INDSIG 

Subroutine to read individual significance test cards and options 



Called from 
WATCEN 

SUBEX? 

SUBEX2 

SUBEX2 

BATHED 

SUBEX2 

SUBEXU 



G3 
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Table VI 



Array 



VARHAM(M^L,MI^F) 
KSTART(MnT) 
TNAME(MNC,MIIT) 
QiNAME(MNCQ) 
X(MNSIGC*80) 
lIYPNAM(MIiL) 
ERRiW.(MIIL) 
OTHMM(MNL) 
LEVSUM( MTJL ,MI^F ) :^LEVSUB 
ILHS(MN}IE) 



Arrays with Variable 
(double precision) 



(double precision) 
(double precision) 
(double precision) 



DUM(MriL,MI^L) 

CCnT(MNL,MNL,MNF) 

;ITABLE(MNP) 

ICQll'Vihrn?) 

:iL£V(MIIF) 

CMPVAR(MNF) 

HDNAM(MNF) 

TRAn(MNL,MNL) 

::?APT(r4NF) 

. r7CUM(MriL,MF) 
:<LC0DE(M1]L,MNF) 
;uIAME(MA,MNF) 

:'i>{K:F) 

IK(MI^F) 

icim?) 
i:j(mnf) 

"!AFCD(MNF,Mr^T) 

ickll(mi^f) 

VA?(MriL) 
: VlLC(MriCD) 

i?L{m?) 

?T.rART(MNT) 

F:'fr(iiFNri') 
'.TrT(M:iL) 

/ r( 3,MII3V,MIICELL) 
::)LnT(r€ICELL,MNF-2) 
, IHM(M^ICELL) 

A(::oRD,:iORD) 

?(:iORD,NORD) 
C(NORD,NORD) 

:;(:;oRD,NOPr)) 

E(NORD,NORD) 



Subscripts 

Subroutines with Fixed Dim ensions 

MAIN,SUBEX1,SUBEXU 

MAIN ,SUBEX1 ,SUBEX3 ,SUBEXU 

MAIN,SUBEXU 

MAIN,SUBEXU 

MAIN,SUBEXl4 

MAIW,SUBEX2,SUBEX^ 

MAIN,SUBEX2,SUBEXU 

MAIN,3UBEX2,SUBEXU 

MAIN,SUBEX1,SUBEXU 

MAIN,SUBEXi^ 

SUBEXl 

SUBEXl 

SUBEXl 

SUBEXl 

SUBEXl 

SUBEXl 

SUBEXl 

SUBEXl 

SUBEXl 

SUBEXl 

SUBEXl 

SUBEXl 

SUBEXl 

SUBEXl 

SUBEXl 

SUBEXl 

SUBEXl 

SUBEXl 

SUBEXl 

SUBEXl 

SUBEXl 

SUBEXl 

SUBEXl 

SUBEXl 

SUBEXl 

SUBEXl 

SUBEXl 

SUBEX2 
SUBEX2 
SUBEX2 
SUBEX2 
SUBEX2 
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Table VI (Cont'd) 
Arrays with Variable Subscripts 

Subroutines with Fixed Dimensions 



Vl(NORD) SUBEX2 

V2(N0RD) SUBEX2 

V3(N0RD) SUBEX2 

V!»(NORD) SUBEX? 

V5(N0RD) ' SUBEX2 

AHYP(M?1CQ) also AHYPJNORD) SUBEX2 

DVEC(MMT) SUBEX3 

E(MAXD) SUBEX3 

IBEG(MKT) SUBEX3 

Nl(MNT) SUBEX3 

H2(t4NT) SUBEX3 

Kl(MNT) SUBEX3 

K2(MNT) SUBEX3 

A(MNL,MnL) SUBEX3 

B{MNL,MNL) SUEEX3 

, C(MNL,MNL) SUBEX3 

D(MNL,MNL) SUBEX3 
NROW(I^T) • SUBEX3 

V(MNL) SUBEX3 

LHE(MNHE,5) SUBEXU 

A(MNL,MNL) SUBEXU 

B(MNL,MNL) SUBEXU 

C(MNL,2*MNL) SUBEXU 

X(MNSIGC*80) SUBEXU 

'"NAME(MNC) GUBEXU 

YD(MNL) SUBEXU 

IROWC(MNHE) SUBEXU 

D(t4NL,Ml4L) SUBEXU 

AINTG(MNL) SUBEXlj 

BDY(MNL-30) SUBEXU 

TOTNAM(MNL) (double precision) GUBEXU 
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Table VII 

Limensioning Constants 

Mr^{L=:]ORD= Maximum number of levels for a factor and one more than the 
maximum number of continuous variables 
Originates in MAi:;, SUBEXl , SUBEX2, SUBEX3, and SUBEXli 
Standard version=^iO 

M:JF= Maximum number of factors (including continv.ous variables) plus 1 
Originates in MAUI and SUBEXl 
Standard version=8 



D=M::T=NM= Maxlm^om number of entries in the data vector; maximum number 
of terms on the significance test card; or maximum number of 
matrices in core at one time (in SUBEX3) 
Originates in MAIN, SUBEXl, and SUBEX3 
Standard veV'sion=258=2^*MNF+2 

M:jC= Maximum number dt words in each term name on the significance test card 
Originates in MAIN, SUBEXl and SUBEXl 
Standard version=3=12 characters 

N!::CELL=- Maximum number of cells for calculation of means and sigmas 
within cells 
Originates in SUBEXl 
Standard versicn=5l6=^*'^*(N-^'F+l )+U 

MAXD= Maxim-ar. di^.ension of a vector in SCP 
Originates in SUBEX3 

Standard ^-ersion=20 ,000 > greater than U^MNL^*2 
(It is importa^ .hat this be is large as possible without 
■ ^ aating the size of storage. It limits the size of the 
san-of-cross-products matrix, '//hen the SCP matrix is not 
large enough to hold the entire set of parameters, two or 
more parses are made in a loop which is the slowest part of 
the program. Making two or more passes through this loop 
.greatly increases the cost of running the program.) 

yjC:i= Maximum number of words in QjlAE^. (a significance test) 
Originates in SUBEXU 
, Standard version=HO=l6o characters 
(Through an error, the printout will only l%t NORD words 
of the test statement.) 

M:ISI';C= Maximum number of significance tost cards 
Originates in SUBEX^* 
Standard version=9 
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Irrensioning Constants 

;FNf:= .'laxinuir. na^ber of words in the variable format statement 
Originates in SUBEXi 
/t^mdard version=*200=10 .^ards 

N!:;!i= Maximum, nunber of items in each of tne LKE : ists 
^rifanates in MAIN anJ i:VhEXl 
'tanuard ver3ion=<?0 

'-CiA= Maxir.ur. number of characters in a factor name (the alphabetic 
iescription on the contrast card) 
Originates in SUBEXi 
L'tandard version=20 

M:irV= Maximum n^amber of variables to be printed across a page for 
cell means and standard deviations 
Ori/^inates in SUBEXI 
v'tanaari version=8 
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Table VIII 
Altei'able Constants 
Unnamed common in MAIN and standard version designations 



KIN 


System input unit = 5 


KOUT 


System printed output unit = 6 


KPNCH = 


System punched output unit = 7 


ISAVE = 


Temporary storage unit = 19 


IThr^IP = 


Temporary storage unit = 18 


ICTOR = 


Temporary storage unit = 17 


IKEEP = 


Temporary storage unit = l6 


I BAG 


Temporary storage unit = 15 


ISIGT = 


Temporary storage unit - 1^ 


lOTFIL = 


Temporary storage unit = 15 > 


ISCPT = 


Temporary storage unit = 12 


ISWT = 


Temporary storage unit = 11 


inFILE = 


Data input file if not KIN 


u:jder = 


Smallest number machine will handle = l.OE-70 


OVER 


Largest number machine will handle = 1*0E70 


CNEARO = 


Constant for controlling rounding error = .000001 


INFO 


Array for controlling analytic sequence 
INFO(35) to INFO (50) unused on June 1, 1972 
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Table IX 
Constants in INFO 



INFO(l) = 1 If analysis of variance 

= 2 If correlation 

= 3 If facto I' analysis 

= 5 It interaction tables 

= 6 If vai-iance reduction 
INFJ(2) = 1 If error sum of squares is residual 

= 2 If error sum of squares is generated from estimates and 
design parameters 
INF0(3) - ANOVA: Degrees of freedom for hypothesis 

CORREL: Number of variables in hypothesis set 
INFO(U) = Order of error sum of squares matri"' 
INF0(5) = A^OVA: Number of variables 

CORREL: Number of variables in error set 
INF0(6) = Degrees of freedom for error 
INFO(T) Controls battery reduction 
INF0(7) = 1 Efroymson's stepwise procedure 

= 3 Step-up procedure i 

= 2 Wherry-Doolitt-le procedure 
I?IF0(8) = Probability, level for adding and dropping variables in 

battery reduction 
INF0(9) and IIIFO(lO) Control standard analysis print-out 
inP0(9) = 0 Print minimal sets 

= 1 Print intermediate sets 

= 2 Print everything 

= 3 Print INFO(lO) output only 
INFOdO) For ANOVA 

= 1 Px^int orthogonalized raw estimates 

= 2 Print disc. fn. coef, for standard scores 

= 3 Print disc, fn. coef. for raw scores 

= h Print transf, matrix for canonical contrasts 
INFO (10) For CORREL 

= 1 Print cetnonical rav score regression weights 

= 2 Print canonical raw score regression weights to obtain unit 
sur.i of squares regressed scores 

= 3 St. score and r^ijw score weights for repressing hypothesis 
variables onto error^variables 
INFO(lO) For factor analysis 

= 0 Print minimal set 

= 1 Print everything 
INFO (10) For VARED 

- 0 Print data for hypothesis variables only 

= 1 Print data about design parameters, contrasts and 
correlations among hypothesis sums of squares 
INFO(ll) through INFOdM Contro] rotation procedures 



y:v-;^_> Nuriber of fcictcrs ir. factor analysis 

^ Criterion for ielern^ini.^g number of significant ca.nonicaI 
variates ^n Mi^ijVA or car^onical correlation rotation problems 
::;F■"^'.. ) = l if Quartinax 
= ? If Yarirr.ax 
= 5 I f Equainax 
= I f Pronax 
= If midtiple f.roup.-^ 
= t If orthogonal centx'oids 
= 7 If orthogonal bounds 
::;rX- = l if factor amvlysls or direct rotation of canonical variates 

^ = If indirect rotation of canonical variates 
I.;pT''l-0 = 1 If ta^onouy of variables 

= ? If taxonomy of groups 
I.\T'Kl3) - .^jr.Der of covariates 
I:;F0(1(0 = 1 Classic^jtl covariance adjustments 

= ? Covariance reiuclion oy error regression weights 
i:iFC(l7) = ! -jr^ber of extension variables 
i:,'FJ'lO = 1 If a test with dimension reduction 

= If a test with aimension reduction and for same 
= ^ If a lest without dimension reduction but for same 
I:JF«^hiq) = :r^-ber of i^mensions which have been set up for redaction 
-^{vQi-'-)) = ::u.T.b'r of observations minus number of parameters svrept out 
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riFOfjl) = : If raw estimates only for AHOVA 
I:;F"<:i; = : rnncipal component.' of dispersion 
= . Principal components of correlation 
= _^ f'rinciDal fact^.r analysis 
= i Image analysis 
= \ Canonical factor analysis 

- ( *!axi urn liKelihAjd 
= 7 /Cpna 
= ?*.inres 

= J comjaunalitif>S' r^-q aired 
= I It' >^onmuna: ities start with S,M,C, 

- / Tf communal I ties supplied 

= Nomb'^r of factors to be remioved and tested = IN?'0(11) 

Cigni f icance ^^st n procedure 
= 0 If none 
= 1 If Pao*s 
= r If Fippe*3 

]p z'^Y) = Pow^^r ^f •! to b^- abol as a convergence criterion (multiplied by 10) 
T ^;,iinber of indivi laa. significance test cards 

:,im.ber of .-^allr to 1111310 (i.e., number of the current hypothesis test) 
rotations are multiple groups, orthogonal centroids or 
r^h'^^r ^r.al b^-^unds 

y^, ( .^^] :iniJ tne n^imber ^ f inu/.d-.al significance test cards read in IV.ICIG 

::V>' ! = >^urren^- VAV^^i ^) , 

f i ^ \ - ii:y' ' H( f.r rr^*vioas toct 

-;i.T ( = [ K- r calls t .j Wi.ks from DINaVKP 

inFn{^,s) - 1 7o prevent copying of data file for reanalysis 
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Malhematicai ^loles 

/tioor rKt«'o art' r. ^ier>.^ribe r.ne mathematical calculations involved 
.^ . ^Iviiv liiicar r.-^iel. Ac the steps are described, the sub/outines 
w>:-i t':o;c.it^. Vner^o 3tcps will be indicated In capital letters within 
: U*ent h»*S«"^o , 

l^io prr-r^L-. is written as a full linear model. That is, as the vector 
f ->a.wrenent3 r?r each observation are read into the computer, dummy 
r;r u-.'lerj. (or iesi.-n parameters or fixed variables) are added to the 
/•_ct-r according lo the requirements of the linear model (DATJEC) . The 
!t*w lata vector Ir.ciulin.^ dummy parameters and measurement variables is 
\ :n tap^- >^r iisk, 

.;....n all .,boervutijnr>hive been read in and data vectors hav^ beer 
•.V./* r^c^ea, \ jam-o*-r>q.iJire3 and cross-products matrix is then constructed 
,\ :r . rujv ^ ) 'uvvTO 'iate analyses which have a large number of 

>. ..-j'-^tr::, t:.- .*> jn-^. f-squares and cross-products matrix is not necessarily 
^ r- i <:.'-ire./ in th- machine. The prof:rain generates and stores as large 
• ^:r.: r-t^rix as p.jssibie, then rewinds itself and generates 
' r^r*: ^r. f ^h«- matrix. Unfortunately when it is necessary 

. . ^ !-.^q..ares and .:roso-products matrices in portions, the 

- 'i r.r. ir.crean^-^. .-,reatly. (lortanately this is a very rare 
- .>'r'':i'^^;. ) Ih- f/ir,-uf-squareo mat t ix is partitioned into submatrices 
V' : w ,rv v^mri^t^^Iv ir-termiiicti by the effects listed on the model card. 
'M^ -x^rix ir j^or^'i m iisi'^ as parMtioncd blocks and elements within 
\. . ratri^-r t r.an elerr.ents within a matrix. 

'.rr:**'- ''i:^- * -.>q larrs '-'i^rix i/- then reduced to the sots -f 
•.-j/irr-r. ani cr-^su pr-d'iCts for th^ various effects as required by the 



-noiel 'iVEEP). This is done by sweeping out the effects listed in the 
model in the order in which they appear ^ from leTt to right) on the 
significance list card, A statement of the orier appears on the print- 
v^ut . The sweeping process is done m st rai^htfoi-ward Gaussian reduction, 
bit instead of reducing the matrix ; -i iiaf^ena] i:. element at a time, the 
process reduces the matr: a part it i red i;lock at a time. That is^ the 
Gaussian sweep is generalized fron an element at a time to a submatrix at 
a ^ime. Otherwise, the process ised is that described by Bock--{ 19^3 ) . 

Subsequent to the sweeping prcc^Gs , ttie" ckppropriate sums-of-squares 
and cross-products matrices for each linear analysis are assembled (TUBEXU) 
and passed on to :he appropriate analytic section. 

AHOVA and CORREL Solutions 

Analysis of variance and correlation problens stex from th*^ determinantal 
equation 

|A-XBi=0 (1) 

where A is a matrix of the sm\ of squares and cross products for hypothesis 

an I B is, for A:inVA, the sum of squares and cross products for errors or, 

for CORBEL, the su.n of squares and cross products for "totals" { isinfr the 

convention that the s\m of squares for errors plus the sum of squares for 

hypothesis is tho sum of squares for totals regardless of what other 

parameters are in the model for the data), 

7h<=> solution of this equation (ANAL) is accomplished as follows. 

1, The matrix B is scaled to correlations and A scaled 
compatibly. 
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■I-he scaled B is decomposed (DECOMP) int.o a product of 
triangular matrices. This triangular matrix is then - 
inverted and pre- and postmultiplied onto the scaled A. 
This results in reducing the original determinantal 
equation into a single matrix eigenvalue problem. 

i'he eigenvalues are calculated and eigenvectors determined (HOW). 

-rius deooription is the essential flow of the calculations. However, 

there are many side statistics calculated to provide expository material for 

the analyses. Table X lists the side statistics and gives their formulations. 

In Table X the symbols have the following meanings. 

Z-the hypothesis variables: in AITOVA the duinray parameters, dummy 
variables design parameters, fixed variables or independent 
variables; in correlation the independent variables or predictors. 
T'hey are assumed to have mean zero. 

X-the error variables: in ANOVA the observed variable., the 
dependent variables, the continuous variables: m correlation 
the criteria or dependent variables. They are assumed to have 
mean zero. 

'^en '-.e determinantal equation (l) is written for correlation problems 
asinc ': mi X, it can be written in matrix notation as 

Kx-:/; (z'x)"^ (z'x) - >(x'X)i = o. (2) 

The v.l-.ec r ^ are the squared canonical correlations. When the equation 
-.J wri^-ton f'-r analysis of variance, the matrix B is 

r = iX'X) - (X'Z) (Z'Z)"' (Z'X) 
„.r.i-h 'XT .v-.nvenience of writing wl.l stir be notated X'X even though 

:z -i rr.3iiual error or wilnin-cel is matrix. The values of X here 
c--::r.'- v-i^ i--3 wf s'lms of squ-irti- 

;' A of equati. r. 1 

i., the ii'ipor.a: mtrix of the square rocts of thp Jiifrnal of 
'/17/ xn correiat i n problems this is divided by a constant 
to Kivf^ the standard deviations of the Z variables. 
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0 is the diagonal matrix of the square roots of the diagonal 
^ of X'X In correlations probU-tas this is dividf'd by a 
constant to give the standard deviations of the ^. In 
analysis of variance this is divided by aconstant to give 
the standard errors of the X. (The ir..itrix X X i.. not 
the same in both solutions.) 

-1 -1 

Q is the triangular decomposition of n^'X'Xo^ 

T is the matrix of eigenvectors of 

iQ'-^^X'2)(Z'Zr^ (Z'X) o^-V^ - Xll = 0 

DFH degrees of freedom for hypothesis 
DFE degrees of freedom for error 
The final matrix equation of the solution can be written 

^-1/2^,Q,-1 -l(X'Z)(Z"Z)-^7.'X)n-^Q-^T'X-^^^ = L 

A 

The various parts of this fonr.uiation are entered in TaMe X .nd cross- 
referenced to Table I to help the reader recognize the printed output if he 
is not familiar with the style of nomenclature used here. 

Factor An al ysis Solutions 

The factor analytic solutions avaiJable in VARAN are common ones and 
the nomenclature used is standard. The subroutines for doing the maximum 
likelihood. ALPHA and MINRES procedures, were written by Michael Browne 
and were merely adapted for use here. These routines print out information 
which i. not generally used but was retained nevertheless. 



Interaction Tables 

The mathematics of calculation of interaction tables is simple enough. 
The parameters for nesting within the cells of the desired interaction are 
generated, and the estimates for this desi.,n ar^- ^-alculated, viz: WAP 
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] r-odu.^es the estimates for the oellb cf tho AP face of -i faoV^rini design. 
pV'^I" < 'near: r.uit nc-^ c^* n":] xdv: ii\ lUlkl-''. »\oU-»;<',, 
Th^ rioronym PJTAPL r^or.-J'ides /OO ^\ - "ilaM onr, vnen e-tir'atc3 
» Hve been 'jV* ai red , th^n prirtj them o^^t, 

V'lriari'^e Peduct ion Analys : r 

/arlance reiuc^ion i.:; i pr'^eriure for evalir^tlng the 

f'l'-^rt of i!i)alar*'?^ in -ui A!i '^'*'- te?^p:n 'lU'. to un* i lal cell oiz'^rs. i^aslcally, 
\\v nypothesis sums of , q iar' ; t v *-he piven tes* \z generated both before 
ani after taking account of t^u? inequality cf evil size^. • (TVie p;rand mean 
is not removed from the data.) The "before*' sun.^ of squar'^s is assun'^d 
to be IQO per cent of the data available. When the "after" suris c>*' -^^quares 
is larger than the "before" a nogatWe less is indicate'!. 

Dimension Reduction Analysaa 

Tnis feature is available throu/?h the significance test card and is 
indicated by a minus sign (-). In use it pertains to a feature of multi- 
variat*^ analysir of variance and canonical correlation which has no uni- 
variate a'^'^^'/F. 

To explain its u'^**, :'ipp.;se w^* are ^?iver two independent ^' tatist i <\al 
nypotheses, both involving the san^ orr<'>r variables, such as an in*-eraction 
mi a main effects test in a MANOVA lesign )f p variables. Tupp^se the 
interaction test, AP-V, has r < p sir.ni Ti'^ant dimensions (or roots), 
rh'^'S^'' r linens ions of th^-^ p dinensi' >n-. ^'f error have the sano 
IV as r variabl^^'^ i ri r 'inivarlat<» an'ilyses and suf^F/^st that 

* f^ain ^^i'Veotr* th'^-s r /ariai:er j • ♦ ^ * x'l- ui'^-i *^U'*Mi*'r • ^iow- 
"V-r, tf)<^» r^maininr p-r ]im'^nr>i' jru*. of fivx' ar*- !i')t int.^»rac^ and c.',ri 
le r'^as'-nably teft'>i "V*^r th* main f^f f^'* -t A. 
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The test of the main effect of A on the p-r dimensions of error is 
reasonable and can be effected by obtaining p-r linear combinations of the 
error variables which are independent of the significant interaction dis- 
criminant variable:; and analyzing these. 

The method of dimension reduction removes the significant interaction 
discriminant variables from the sums of squares for the test A=V, This 
forces those same linear combinations of the error variables to have zero 
roots in the determinantal equation |A-XB| = 0, thus removing them from 
the analysis. The main effects test A=V is then calculated, and the 
multivariate probability tests are compiled as if there were only p-r 
variables in the test. The univariate F ratios are not altered. 

The notation on the significance test card to produce this test of A 

is 

AIIOVA : W=0 , A=V-AB=V , 
with an individual significance test card used to indicate the probability 
level in the test of AB=V which determines both significance and the 
number of canonical variates to be removed frcm A=V. 

In addition to this application to analysis of variance, dimension 
reduction can also be applied to canonical correlation, 

COPREL : W=0 , V1=V2-V3=V:' . 
-mi -13 'i type Of analysis of covariance, 
i\:;nyA :W=0 , A=V1-V?--V1 . 
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Rotation Procedures for Canonical Variates and Factors 

Eight rotation procediores are available for multivariate solutions to 
linear models. Four of these are familiar techniques of factor analysis: 
Quart imax, Varimax, Equamax and Promax. These four techniques apply 
mathematical criteria to the solutions and rotate the correlations between 
the canonical variates and the variables to a solution that is presumably 
more comprehensible to the researcher. 

The remaining four rotation procedures apply external criteria to the 
rotational procedure. In multiple groups, orthogonal centroids and orthog- 
onal bounds, a set of regression weights are entered and used to determine 
the correlations between the rotated canonical variates and the variables. 
This makes all these procedures Procrustian in form; however, they are not 
explicitly Procrustian in that they use regression weights, not target co- 
efficif^nt matrices. The multiple groups procedure obtains an oblique 
solution. The orthogonal centroids scheme is similar but uses the first set 
of weights optimally, then uses the i^econd set of weights optimally but 
orthogonal to the first set. Successive weights are fitted oj^timally but 
orthogonal to all the previously fitted weights. 

Orthogonal bounds is Procrustian in form also. First a multiple group 
ODlique solution is found then an orthogonal envelope is found which has a 
least squares fit about the multiple group solution. The correlations 
between the canonical variates and the variables are rotated to the envelope. 

Th*^ orthogonal Procrustes scheme is due to Norman Cliff and is 
referenced in the bibliography. The calculation scheme is due to Roger 
r-nmaJ. Hpre the weights supplied are the "target" matrix of coefficients 
and not regr^-ssion weights. 
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CORReiATIUriS bETwEEN THE ERKUK VAKlAbLES AND THE 
KUTaTEO canonical VAKlATtS 

CANOivilCAL VARIATE 
VARIABLE 12 3 



ERROR 1 
ERROR 2 
ERROR 3 



0*9916 
-0*061b 
-O.lZbO 



-p*0394 -O, 1149 
0*96bt> 0*24^3 
0.2464 0*9611 



TRANSFORj'iATIUi r.ATRIX 

Tka.^6. VECTOR 

CAiM* VAR. 12 3 

0.6B0b -0.4192 0,6007 

-u,2v9^ 0,3892 0,7b0b 

0 • 66 y 0.6907 -0 • 2 V b6 

i>UK-Ut— SQUARES c.aTRIX OF RUTaTELi HYPOTHESIS CANONICAL VaRIATES 

CAjnjUNICaL VARIaTE i 

12 3 ^ 

1 0.1151 -0.0747 0.0669 

2 -0.0747 0.07o4 -0.01b7 

3 0.0669 -0.0137 0.1333 



'.ieiOHTS FOR REGKtSSlHG ERROk VAklAr^LbS OivTO ROTATED ERROR CAi^DMCAL VARlATFb 



CAimOimICAL VARIATb 
VARIABLE 1 2 



ERROR 1 
ERROR 2 
ERROR 3 



0.0063 
0.0017 
0.0047 



0. 0002 
0.0r>32 
-0.0113 



3 

0. 0008 
-0.0139 
0.0467 



TEST WITH fiO ERROR }n*GT DONE 

A = 0 

TEST WITH NO ERROR NOT DONE 
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V/-\RAAi: LINEAR MODEL VAIUaimCE ANALYSIS 
FIRST edition: JUNE 1, 1972 

fducatiomal testing sprvice 
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